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Abstract: We provide a rigorous and asymptotically exact expression of the mutual 
information of translationally invariant free fermionic lattice systems in a Gibbs state. 
In order to arrive at this result, we introduce a novel framework for computing deter- 
minants of Toplitz operators with smooth symbols, and for treating Toplitz matrices 
with system size dependent entries. The asymptotically exact mutual information for 
a partition of the one-dimensional lattice satisfies an area law, with a prefactor which 
we compute explicitly. As examples, we discuss the fermionic XX model in one di- 
mension and free fermionic models on the torus in higher dimensions in detail. Special 
emphasis is put onto the discussion of the temperature dependence of the mutual in- 
formation, scaling like the logarithm of the inverse temperature, hence confirming an 
expression suggested by conformal field theory. We also comment on the applicability 
of the formalism to treat open systems driven by quantum noise. In the appendix, we 
derive useful bounds to the mutual information in terms of purities. Finally, we provide 
a detailed error analysis for finite system sizes. This analysis is valuable in its own right 
for the abstract theory of Toplitz determinants. 
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1. Introduction 

How do the correlations of natural quantum states of many-body systems behave? If 
"natural" is taken to mean "generic" in the sense of a random pure state drawn from 
the Haar measure, then the answer to this question is: Subsystems will almost surely be 
very nearly maximally correlated with their complementary subsystems. However, this 
is not the situation that one is usually interested in in many -body and condensed-matter 
physics. Ground states of local Hamiltonians typically exhibit far less entanglement 
than that suggested by the previous argument. Indeed, there is a large body of evidence 
|[TU33lL suggesting that ground states of gapped quantum many-body systems on lattices 
satisfy an area law 1 1 ] for the entanglement entropy. In other words, given the pure state 
of a lattice system, if one distinguishes a certain (connected) region of the lattice, then 
the von-Neumann (or Renyi) entropy of the reduced state associated with this region 
does not grow, as one might expect, like the number of degrees of freedom of this region 
(its "volume"). Instead, it scales like the number of degrees of freedom on the boundary, 
hence as its "area". Area laws have been shown for all gapped one-dimensional models 
with local interactions (24l [30] [31] . Extensions to classes of higher-dimensional lattice 
models have also been obtained (21 [12J |29) . These results and others therefore suggest 
that ground states of quantum many-body systems are much less entangled than they 
could be. These measures of entanglement - refined quantities revealing much more de- 
tailed information about the structure of correlations than more conventional correlation 
functions - in a way inherit the decay of correlations. This deep insight is also at the ba- 
sis of classical efficient simulations of quantum many -body systems: The formulation 
and analysis of many-body systems in terms of matrix-product and other tensor network 
states have put this intuition on a solid theoretical footing (T][34][35l. In numerical sim- 
ulations of many-body systems, the importance of such approximations can hardly be 
overestimated. As a consequence, the last decade has seen an enormous amount of in- 
terest in the study of entanglement in quantum many-body systems in the condensed 
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matter context and in issues of numerical simulation fTl l34l[36lL entanglement spectra 
l36lfiTlL and their relationship to quantifiers of topological order 138 -411. 

Can this analysis be extended to thermal states? A moment of thought reveals that the 
entanglement entropy should in fact satisfy a volume law for thermal states. However, it 
turns out that the entanglement entropy is not the right quantity to grasp the correlations 
of a mixed quantum state. There exist many such measures, but the most natural and 
the most frequently adopted one is the mutual information. For a subset A of sites of 
the lattice and its complement B, it is defined as 

I(A:B) = S(p A ) + S(p B )-S(p), (1) 

where pA and pB are the reduced states of p with respect to A and B, respectively. 
For pure states, the mutual information reduces to twice the entanglement entropy. One 
might then hope that this quantity fulfills an area law for thermal states of local hamil- 
tonians on a lattice. This turns out to be true in fact for any fixed temperature I42ll43l . 
Specifically, for the Gibbs states of a Hamiltonians with local interactions on a lattice, 
and with bounded operator norm \\h\\, one finds that 

I(A:B)<20\\h\\\dA\, (2) 

where dA denotes the boundary area of the region labelled A. This statement is in fact 
surprisingly simple to prove. It following rather directly from the extremality of Gibbs 
states with respect to the free energy. 

The strength of this result - that it is completely general - also constitutes its weak- 
ness; Eq. ^ does not say anything about the correlation behaviour of thermal states of 
specific models. Given that the bound is linearly divergent in the inverse temperature, 
the law becomes less and less tight in the limit of small temperatures. In particular, at 
zero temperature, we know that there can be logarithmic corrections to the area law 
for critical systems. Hence, it would be desirable to have bounds which depend more 
explicitly on the properties of the system at hand. What is more, asymptotically exact 
results for important classes of "laboratory" models seem important. 

In this work, we present asymptotically exact results for important classes of "labo- 
ratory" models. We introduce a framework capable of rigorously computing the mutual 
information in free fermionic lattice models. This framework is based on new approx- 
imation techniques as well as on proof tools for dealing with Toplitz matrices with 
smooth symbols. We separate the presentation in a discussion of the one-dimensional 
case and one of the higher-dimensional situation of a cubic lattice with the topology of 
a torus. We find rigorous area laws for the mutual information that are meaningful in 
the limit of small temperatures, and converge to the known results for the entanglement 
entropy for ground states ['10 , 26 ]. This feature is even present in higher-dimensional 
fermionic lattice models fT5ll2Tl - [23ll25ll . For low temperatures, we find a logarithmic 
divergence of the mutual information in the inverse temperature, in fact confirming 
predictions of conformal field theory. This scaling suggests that in order to detect crit- 
icality, extremely low temperatures are needed. This is an exponential improvement in 
the inverse temperature to the above bound. 

Mathematically, we discuss new methods of approximation, allowing to use Toplitz 
matrix techniques to studying arbitrary subsystems of translationally invariant fermionic 
models. What is more, we introduce a novel second order Szego theorem, simplifying 
Widom's formula for smooth symbols. This approach allows to use the trace formula 
without encountering an infinite sum leading to arbitrarily highly oscillatory terms. 
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Hence, apart from the application in the study of quantum many-body physics, we ex- 
pect our results to be a highly useful tool in the abstract theory of Toplitz determinants. 
The paper is structured as follows: 

- In Section [2] we introduce Renyi entropy mutual informations for free fermionic 
models, discuss fermionic covariance matrices and present new, easily computable 
and practical upper and lower bounds to mutual informations. 

- In Section [3] we go on to formulate the problem, and we provide a synopsis of the 
mathematical argument. A general proof strategy for computing the mutual infor- 
mation is also provided. 

- Section [4J] contains the main theorem - the asymptotically exact expression for the 
mutual information - and the core of the proof. Since the argument is general enough 
to capture expressions deriving from all Renyi entropies, the complete knowledge of 
the spectrum of reduced states is also obtained in this fashion. This section reports 
also the main technical progress when it comes to dealing with Toplitz matrices with 
smooth symbols. 

- Section[5]is dedicated to a thorough discussion of the temperature dependence of the 
area law, with an emphasis on the the very low temperatures behavior. 

- We then turn to higher-dimensional free fermionic models on the torus in Section [6] 
which we can treat with similar methods. 

- Finally, we present an outlook in Section [7] comparing the findings with predictions 
of conformal field theory, and discussing implications to the study of entanglement 
spectra l36lfiTl . We also consider the implications for noise-driven, open fermionic 
quantum many-body systems, as have recently been discussed also in the context of 
open Majorana wires and notions of noise-driven criticality l44l[48ll . 

2. Free fermionic models 

In this work, we focus on isotropic translationally invariant fermionic models. For now, 
we assume the most general form for the couplings. Later, we will concentrate on the 
fermionic variant of the XX spin model (under the Jordan- Wigner transformation), 
which is a particularly important special case. In this section, we introduce the class 
of models discussed here. 



2.7. Hamiltonians. We consider free fermionic models on cubic lattices (7Ln) d for 
some even TV G N as TV becomes large. Bi-sected geometries on the torus can be 
related to the situation of having D = 1. For simplicity of notation, but without loss of 
generality, we will therefore present this one-dimensional case in the main part of this 
work, while discussing the higher-dimensional situation in Section [6] The Hamiltonian 
takes the general form 

H: = E fjVjfifk, (3) 

with V = V T G TR, NxN being a circulant matrix, referred to as the Hamiltonian 
matrix. We assume the coupling to be an specified, finite-ranged interaction of arbitrary 
interaction length r G N. The couplings are defined by a sequence of numbers (vk)keN 
with the property that = for k > r. Then take 



dj := v\j+i\ : j = —TV/2 + 1, . . . , TV/2, 



(4) 
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B 

Fig. 1. Geometry of the considered situation for one-dimensional system, consisting of N degrees of free- 
dom, coupled with a local Hamiltonian equipped with periodic boundary conditions. The distinguished region 
embodying sites {0, . . . , L — 1} is referred to as A, its complement is called B. 



having a natural reflection symmetry, d-kmodN = dk, and consider the Hamiltonian 
matrix with entries 

Vij = d(i-j) modiV (5) 

for i,j G Sat. This situation is depicted in Fig.[T] The most important special case is 
constituted by the fermionic variant of the XX model, which is originally a model for 
spin- 1/2 systems. Here, we directly consider its fermionic instance, by virtue of the 
Jordan Wigner transformation. In this language, 

vi = a, v 2 = b, v k = (6) 

for a, b G R and 2 < k G N, so that V is a circulant matrix with a on the main 
diagonal and b on the first off-diagonal. All these models are integrable and exactly 
solvable. Thermal states are defined entirely by the collection of second moments of 
fermionic operators. The arbitrary finite range of interactions has in the main text been 
chosen for simplicity of notation only, and exponentially decaying interactions can eas- 
ily be accommodated as well. A generalisation to exponentially decaying interactions 
is possible. Since V is circulant, we find the spectrum {ek} of V to be 

1 J2 e k e^ k / N = d k . (7) 
ke% N 

For the well-studied fermionic variant of the XX model, see Eq. ([6]), one gets the famil- 
iar expression 

/27Tife\ 

£ /c = a + 26cos( -jj-y (8) 



2.2. Majorana fermions and covariance matrices. It is convenient to define fermionic 
covariance matrices in terms of Majorana fermions: 

f. i ft f . _ f\ 

n := r l+N := ^p-, i = 0, . . . , N - 1, (9) 

similar to the canonical coordinates for bosonic operators. These Majorana fermions are 
Hermitian, traceless, and form a Clifford algebra. We are interested in thermal states of 
the quadratic Hamiltonians, which are particular instances of (quasi)-free or Gaussian 
fermionic states. Such states are completely specified by their fermionic covariance 
matrix |25l|44l|49l|5l7GR 2Ar><2Ar with entries 

7i,j = itT (p[ r i, r j]) > ( 10 ) 

i,j = 0, . . . , 2N — 1, where the brackets denote the commutator. 7 is always anti- 
symmetric, 

7 = -7 T , (11) 

and satisfies —-f 2 < 1. 
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2.3. Covariance matrices ofGibbs states and their reductions. Matrix functions of the 
covariance matrix can be computed exactly. This is the case because any covariance 
matrix is unitarily equivalent to a direct sum of 2 x 2 covariance matrices, reflecting a 
situation of entirely uncoupled fermionic modes. This is the fermionic analogue of what 
is often called the Williamson normal form in the bosonic setting. As a consequence, 
one can identify an explicit expression for the covariance matrix of Gibbs states 



e -PH 



(12) 



r tr(e-f> H ) 

at inverse temperature fj > 0. The covariance matrix 7 G R 27Vx2Ar of p is given by 

7-f ° f{h) ) (13) 

where / : R — >• R is defined as 

fix) := 2^^ 1 = - tanh ( ^] . (14) 

We suppress the temperature dependence here: Throughout this work, we will be con- 
cerned with Gibbs states with respect to some /?. We allow for arbitrary temperatures 
and will also later consider the asymptotic limits f3 — > 00 and /3 —> 0. 

Reduced states of Gaussian states are always Gaussian (as can be seen most easily by 
considering their Grassman representation; compare Ref. [49]). The covariance matrix 
of a reduced state is the appropriate principle sub-matrix of the full covariance matrix. 
Applied to Gibbs states, we find that the reduced state of a subset A = {0, . . . , L— 1} c 
{0, . . . , iV — l}isa Gaussian state with covariance matrix 

7U = f_„%. /<A » U l, 05) 



-f(h)U 

where .\a denotes the principal sub-matrix associated with the degrees of freedom of 
subsystem A. 



2.4. Toplitz matrices. Toplitz matrices I5TH531 may be viewed as principal sub-matrices 
of infinite circulant matrices. Such matrices will play a crucial role in this work. Through- 
out, we will encounter families of Toplitz matrices T n G R nxn , whose entries are given 
by 

(T n )i,j = U-j (16) 

for some sequence of reals (tk)kefi- Colloquially speaking, such Toplitz matrices largely 
resemble circulant matrices, with the "upper right and lower left corners" deviating from 
a strict circulant matrix. In fact, most of the theory on Toplitz matrices (Toplitz deter- 
minants) is in one way or the other concerned with the error made when replacing a 
Toplitz matrix by a circulant matrix of the same dimension. Later on we will consider 
sub-matrices of large finite circulant matrices. It is clear that this sequence of Toplitz 
matrices is defined entirely in terms of the sequence of numbers (tokens- This sequence 
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of numbers - and hence the sequence of Toplitz matrices - is most conveniently repre- 
sented in therms of its symbol, defined as the inverse Fourier transform on the torus 
T := {x G CD : \x\ = 1} of (tk)ke%> i- e - me symbol t G L^T^is defined as 

th = h C det{ei9)e ~ lk9 ■ (17) 

The decay of the Fourier coefficients is directly related to the regularity properties of 
the symbol. The summability of the coefficients (i.e., ^ ke ^ \tk\ < oo) is sufficient to 
ensure t G Loo(T). The relationship between the decaying behaviour of the Fourier 
coefficients and the regularity of the symbol will be used very frequently. In fact, given 
our assumptions on the interaction parameters, much better regularity results can be 
derived as we will show and use later on. Mathematically, t G ^oo(T) is equivalent 
to the requirement that the associated Toplitz operator generates a continuous linear 
operator on / 2 (N). Furthermore, note that all spectra can be expressed entirely in terms 
of the symbol. 

The starting point of our analysis is the following observation: covariance matrices 
of subsystems of thermal states of translationally invariant fermionic models are well 
approximated by Toplitz matrices. Since the Hamiltonian matrix V in Eq. ^ is circu- 
lant, the covariance matrix 7 of a Gibbs state p is, for any fj > 0, circulant as well. That 
is to say, 

7i,j = k-jmodN, (18) 

for some real sequence (lk)ke%> suppressing the temperature dependence. Indeed, sub- 
matrices j\a for a region A — {0, . . . , L — 1} of the full covariance matrix 7 (of the 
large but finite system) can be well approximated by the Toplitz matrix 7l G R LxL for 
large TV (see Section [42] ). Th family of Toplitz matrices 7^ can be expressed in terms 
of the symbol 

X = foe (19) 

where e e -Loo (T) is defined as 

e{x) = Y,d k x k . (20) 

Again, via the inverse Fourier transform on the torus, one recovers 

d k = ^~ 1^ d0e(e w )e- [ke , (21) 
^ Jo 

where (dk)ke% is j ust me sequence of numbers that govern the coupling in the Hamil- 
tonian matrix in Eq. ([3]). 



1 As usual, L p (T) denotes the set of equivalence classes of p— integrable functions up to functions that 
vanish almost everywhere on T, with p = 00 denoting the equivalence classes with finite essential supremum. 
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2.5. Entropies offermionic Gaussian states. Given that the states we consider are com- 
pletely described by their covariance matrices, we are able to explicitly compute their 
Renyi, and specifically von-Neumann, entropies explicitly and efficiently. In particular, 
we consider covariance matrices of the form 

where X = X T G R nxn . For fermionic Gaussian states p of n modes with a covariance 
matrix of this form one finds 

Sa(p) = ^— log 2 tr(p") = tr( 5a (X)), (23) 
1 — a 

where the function s a : [— 1, 1] — >> [0, 1] is defined as 

Taking the limit as a —> 1, one recovers the von-Neumann entropy, 

S(p) = tr(s(X)) (25) 

with s : [-1, 1] -> [0, 1] being 

These expressions will be central for our analysis. 



3. Statement of the problem 

J. 7. Computing mutual informations. We will now turn to the main object of our study, 
the mutual information', a measure of the correlations between two non-overlapping 
subsystems. We will consider a one-dimensional lattice with N lattice sites and let 
A = {0, . . . , L — 1} constitute one part and B = {L, . . . , TV — 1} its complement. The 
quantum mutual information between A and B is defined as 

I(A : B) = S( PA ) + S(p B ) ~ S(p AB ), (27) 

where again, p is a state on the entire system, pA (Pb) is obtained by tracing out sub- 
system A (B). The mutual information naturally captures all correlations, quantum as 
well as classical, between A and B. It is a meaningful measure of correlation for mixed 
states, including finite temperature thermal states as a special case. For zero tempera- 
ture, it reduces to (twice the) entanglement entropy. Similarly, expressions of the kind 

I a (A : B) = S a (p A ) + S a (p B ) - S a (p) (28) 

for an a > different from 1 are well defined mutual Renyi entropies. We will identify 
novel formulae for the asymptotic behavior of the mutual information for Gibbs states 
of isotopic translationally invariant free-fermionic models, and present bounds that al- 
low to study the limit of large f3 analytically. The formulae given will be exact in the 
asymptotic limit of large TV. 
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3.2. Structure of the argument for one-dimensional systems. A number of steps will be 
necessary in order to arrive at an asymptotically exact expression of the mutual infor- 
mation. To start with, we need to have a handle on how to make the intuition rigorous 
that we can compute S a (pA) (and S a (ps)) as if it was the reduced state of an infinite 
system. The result has to fit then with the expression for S a (p) for the entire system, in 
a way that only the boundary terms remain. The turns out to be a delicate a faire, and 
requires very different tools than the pure state analysis, where the asymptotically exact 
entanglement entropy is obtained using the Fisher-Hartwig formalism [10]. In the case 
of thermal states, the so-called "double scaling limit" makes sense, where A is taken as 
a constant fraction of the total system size, and the total system is taken to infinity. 

We aim at computing an asymptotically exact approximation of the entropy of a 
subsystem 

S a (p A )=tT(s a (X\ A )), (29) 
where X G x N is given by 

X Kl := x k ^ := 1 X^'^e-^-W. (30) 

3=0 

The proof strategy is as follows: 

- We perform a continuum limit on the full system 

4? Xif ■= 4-i == i fj H^)e-^fd4>, (31) 

and show that the entropy of a subsystem can be computed with an exponentially 
small error in N. 

- When computing entropies of sub-matrices S a (pa), we can in this continuum limit 
invoke the theory of Toplitz matrices, even though the entries of the Toplitz matri- 
ces actually depend (very slightly) on the system size. We can therefore consider 
families of Toplitz matrices with symbol 

\ = foe. (32) 

The same approach is feasible for S a (pB)> 

- We apply trace formulae of Toplitz matrices with smooth symbols. This will allow 
us to compute an asymptotically exact expression of the entropy of both subsystems. 

- We find an expression for the entropy of the total system S a (p), asymptotically exact 
in the limit of large N. 

- The expression for the mutual information found in this way will contain an infinite 
sum, reflecting the infinite number of modes in momentum space in the asymptotic 
limit. As such, this formula can not be evaluated. This obstacle we will overcome 
by a new technique that we introduce, making use of the theory of highly oscilla- 
tory integrands. This technique is expected to be of significant interest also outside 
this context, when analysing properties of families of Toplitz matrices with smooth 
symbols. 

The combination of these steps will allow us to provide an asymptotically exact rigorous 
and computable expression for the mutual information. 
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4. Main result 

In this section, we will state the main results for one-dimensional systems. We will first 
present the main theorem and then continue with the proof of the theorem. This will 
require a number of techniques that we will lay out in later subsection s. The specifically 
important case of the XX model, the temperature dependence as well as the situation of 
free fermions on the torus will be discussed in later sections. 

In accordance with the analysis of the previous section, we will evaluate the mutual 
information by analyzing the determinants of the covariance matrices, as the system 
is taken to the thermodynamics limit. Covariance matrices of free fermionic systems 
in one dimension are Toplitz matrices, so we are free to use the tools available from 
the theory of Toplitz determinants. In particular, it is well known that Toplitz determi- 
nants behave in a very regular manner as the dimension of the matrices are taken to 
infinity. In the case when the symbol is continuous, Szego's strong limit theorem l57l 
gives the precise scaling of the determinant with the dimension of the Toplitz matrix. 
When the symbol has discontinuities, then it is necessary to use the Fisher Hartwig for- 
mula in order to obtain accurate asymptotics. As we will be dealing with thermal states, 
the symbols will always be continuous, even though a discontinuous symbol reflect- 
ing the Fermi surface can be arbitrarily well approximated for low temperatures. This 
creates a quite intriguing situation: We can "approximate" the situation covered by the 
Fisher-Hartwig theorem arbitrarily well with smooth symbols, and can hence "interpo- 
late" between these situations. As mentioned before, the methods developed here are 
expected to be useful also outside the context of quantum many-body systems. 



4.1. Main theorem for one-dimensional systems. We consider the mutual information 
in the large TV limit, where 

\A\=L = 9(N), \B\ = 0(N), (33) 
meaning that both subsystems grow linearly in N. Specifically we assume 

A = {0,l,...,\qN]-l} (34) 

and 

B = {\qN],\qN]+l,...,N-l}. (35) 
The main result for one-dimensional bi-sected systems can be stated as follows. 

Theorem 1 (Asymptotically exact expression for the mutual information). For any 

inverse temperature (3 > 0, and for any a G [0, oo), the mutual information is given by 
the, in N asymptotically exact, expression 

K ' 47r 2 i_ 7r v J_ n A(e lf >) - A(e^) tan ((</>- 0)/2) W 

(36) 

Here, the temperature dependence is only implicit in A = / o s (since / depends on /?), 
while A' denotes the derivative of A, which has to be understood in the following way, 

\'(6) := ^A(e ie ). (37) 

This is an asymptotically exact and simply evaluated expression. The most important 
instance is the one for the von-Neumann mutual information. 
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Proposition 1 (Mutual information as an infinite sum). For any inverse temperature 
P > 0, and for any a G [0, oo), the mutual information is given by the, in N asymptot- 
ically exact, expression 

x sin (k(0 - (j))) {\'{<j>) - A'(0)) + o(l). (38) 
The proof of this statement will require some preparation. 



4.2. Approximation statements. In this subsection, we collect approximation statements 
that are being used in the proof of Theorem [T] We have a more detailed look at the ex- 
pression 

N-l 

X i N i ] = 4-/ = jj E \{e 2 ^' N )e- 2 ^ k - l ^ N (39) 
i=o 

as well as at 2 

*M = x k ^ = ^ / A(e^)e-^-^#, (40) 

where the entries of the finite Toplitz matrices are always assumed mod N. We know 
that 

Sa(PA) = tr( Sa (I W U)), = tr(s a (X^\ B )). (41) 

Expression Eq. ( [3T] ) is the continuum limit of Eq. ( [30| ). Therefore, from a physical 

perspective, it seems justified to replace x\^_] by its infinite- system counterpart x^j 
in the calculation of the entropy if the system is large enough. In the next lemma we 
will show that for any fixed temperature, the error made by this replacement decays 
exponentially in TV. 

Lemma 1 (Approximation of entropies of subsystems). For A, B given by Eq. 
and Eq. p5| ) the errors 

e A: p(N) :=\tr(s a (X^\ A ))-tv(s a (X\ A ))\ (42) 

and 

e B A N ) ■= I tr(s a (X^\ B )) - tr(s a (X\ B ))\ (43) 

are exponentially small in N for fixed f3 > 0. More precisely, there exist some N -inde- 
pendent constants ola,p,&b,p > 0, (depending on ft, q and the explicit form of the 
symbol A ), such that 

lim ex-p(Na x ,p)e k A N ) = for X G {A, B} . (44) 

N^oo 

Proof. The proof for the set B is completely analogous to the proof for the set A, so we 
will consider the latter subset A only. We will proceed in four steps: 

1. We use regularity properties of s a to bound the total error by the error in the eigen- 
values originating from the replacement of X^' \a by X\a- 

2. We bound the latter by the square norm of the difference of both matrices. 

3. We write the square norm as an appropriate sum of — x k \. 



12 



H. Bernigau, M. J. Kastoryano, and J. Eisert 



4. Note that — Xk | is the well-known error of approximating an integral by a sum 
using the trapezoidal rule. 

For technical reasons, it is adventageous to reverse the order. So we start with the last 
step and consider step three, step two and step one afterwards. The function <p ^ A(e 1( ^) 
is real-analytic and 2tt— periodic with values in some compact subinterval of (—1,1). 
By the first part of Lemma [2] there exist some constants k > and M > such that 



where 

This implies 



A (e i0 ) < M for all e D, 
D:={xe€: \S(x)\ < k} . 



|A(e^)e ifc ^| <Me |fc|K 
for all k G 2. Therefore, by Lemma [2] we get 



(TV) _ 

For the third step note that 



2Me^ a 2Me La _ 2Me^ a 

- e aN _ i ~ e aN _ J e a7V - 1 



(L-l) 

J=-(£-l) 



(AT) _ 

x j x j 



(45) 
(46) 
(47) 

(48) 
(49) 



which follows directly from counting the number of entries in both Toplitz matrices. 
For the second step let fj,r denote the r-th eigenvalue of and let fj,r N ' denote the 
r-th eigenvalue of X^ . Possibly relabeling these eigenvalues, the Hoffman- Wielandt 
theorem (Lemma|3]i yields 

L 1 o o 



E 

r — 



U (N 7 A) _ ,,(A) 
rr H'r 



< 



xW\a-X\ 



(50) 



Finally, note that the Renyi entropy functions s a are Holder continuous (see Lemma]?]) 
with some appropriate Holder exponent < 7 < 1. In other words, there exists some 
constant C > such that 



\s a (x) -s a (y)\ < C\x-y\ 7 

This gives the estimate 

L-l 

caAN) <^\s a U N ^)-s a U L) ) 

r=0 

< c (j2 - 4 L) 



(51) 



(52) 



\r=0 

< c (j2 U^) - 

\r=0 



7/2 



L l- 7 /2 
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where we have used Holder's inequality in the third line. Combining Eq. (52), ( [50] ), 
(49k and ([48]), the desired result follows for any constant 



(53) 



Moreover, since 7 can be chosen to be 1 for a > 1 and arbitrary closed to 1 for the 
von-Neumann entropy, any rate 



ola,p < «(1 - q) 



(54) 



will satisfy the required condition for a > 1. Repeating these steps for subset B yields, 
for any rate 

<*B,p < 7™2S (55) 

the validity of the lemma. □ 

Lemma 2 (Error estimate (54)). Let h : R — > R be real analytical and 2ir-periodic. 
Then there exists a strip D = Rx(- k)cC with k > such that h can be extended 
to a holomorphic and 2it -periodic bounded function h : D — > CD. 77z£ error for the 
rectangular rule can be estimated by 



< 



4irM 



(56) 



where M denotes a bound for the holomorphic function h on D. 



Note that the constant k will not depend on n, however it will clearly depend on the 
function h. It is natural to ask how k will depend on f3 if M is kept at a fixed value 
when we choose h to be the symbol of the covariance matrix of the thermal state of, 
say, the fermionic instance of the XX model (see Section^. The answer is that k oc 
for large f3. 

To see this fix some a > and define 



M a := sup I |tanh(x) \^(x)\ < | j 



(57) 



Note that M a < 00 if and only if a < 1 since the hyperbolic tangent has a pole at ±i| . 
Define the set 



D a :={ze€\\% [p (| + 6cob(z))) j < ^ } 



(58) 



then obviously h(z) < M a for all z G D a . From ^ (cos(z)) = sin(5R(z)) sinh(— $s(z)) 
follows: 

/ a7r \ 



D a := j 



1 9 (z)| < arsinh 



V26/?y 



(59) 



For fixed a < 1 set k := arsinh and D := D a . Then ft, is bounded by the 

positive real number M a on D. Moreover for any a > 1, ft is necessarily unbounded 
on The asymptopic scaling in /3 follows from 
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Lemma 3 (Hoffman and Wielandt [55 |). Let A,Ee M n be normal matrices and let 
(Ai, A2, • • • , A n ) be the eigenvalues of A in some given order and let (/ii, /i2? • ' ' •> Mn) 
be the eigenvalues of A + E in some order. Then there exists a permutation a G S n 
such that 



(0 



2 • 



(61) 



Lemma 4 (Holder continuity of entropy functions). The (covariance) Renyi entropy 
functions s a and the (covariance) von-Neumann entropy function s are Holder contin- 
uous on [ — 1, 1], i.e., there exist < 7 < 1 and C > such that 



\s a (x) -s a (y)\ < C\x-y\ 7 for all x, y G [ - 1,1]. 

Possible Holder exponents are 

- < 7 < a /or a < 1, 

- < A < 1 /or von-Neumann entropy and 
-0 < A < 1/ora > 1. 



(62) 



Proo/ To start with, note that the first derivative of s a (see Eq. ( |141| ) in the Appendix), 
is uniformly bounded in x and y whenever a > 1, which proves Holder continuity with 
Holder exponent A < 1 for a > 1. For < a < 1, there exist constants Ci, C2 > 
such that 



\s' a (x)\ < (Ci + ^ ((1 + + (1 " • 



This yields for any — 1 < x < y < 1 



\s a (x) - s a (y)\ 



s' a (z)dz 



(63) 



(64) 



< 



= Ci \y - x\ + C 2 (|(1 + y)°> - (1 + x) a \ + |(1 - 2/)" " (1 - *) a |) • 



Therefore, 



|s a (a?) - s a (2/)| 

liniSUp : < C 2 . 



(65) 



Hence, for every x G [0,1], there exists some open neighbourhood U x C [ — 1,1] such 
that 

|s a (x) - s a (2/)| 



12/ — ^1 



<2C 2 



(66) 



for every y £ U x . Then Holder continuity follows from compactness of [ — 1, 1]. This 
proves the validity of the statement for a < 1. For the von-Neumann entropy note 
that by Eq. ( 143 ), s' has logarithmic poles at ±1. Following the same line of reasoning 
again yields Holder continuity of the (covariance) von-Neumann entropy for any Holder 
exponent A < 1. □ 
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4.3. Relating the mutual information to trace functions of Toplitz matrices. Given that 
we are interested not in the entropy of a subsystem, but in the mutual information, we 
will need to consider a second order asymptotic formula for the entropy functionals. 
This is due to the fact that the bulk contribution of the entropy is cancelled out in the 
expression for the mutual information, and we are left only with the contribution orig- 
inating from the boundary. We will therefore require a Szego's strong limit theorem 
which contains explicit expressions for the second order contributions. In order to for- 
mulate the theorem, we will have to introduce some definitions and notation. It will 
be necessary to consider sequences of truncations of Toplitz matrices - in our case 
principal sub-matrices of covariance matrices. We will define the following classes of 
symbols. The Wiener algebra 



W := {n € Loo (T) : 



i r 27r 

hi 



j d9 



< CO 



(67) 



the Besov space 



Bj:=UeL 2 (T):£(|n| + 1) 



± n(e w )e- ine d9 



< OO 



(68) 



the Krein algebra 



(69) 



and the space of piecewise continuous symbols, denoted by PC (T) play an important 
role; see Ref. [53 ] for more details. Consider a continuous symbol A : T — > <D and a 
point x range (A). There exists a (unique up to some constant offset of the form 2nk 
with feel) continuous argument function 



arg(A - a) : [ - 7r, n] ->> R; <j) i-> arg [A(e 1( ^) - a] . 



(70) 



Independent of the choice of offset, the winding number 



wind(A, a) := ^— arg(A — a)(7r) — arg(A — a)(— tt) 

Z7T 



(71) 



is a well-defined integer. Before formulating Szego's limit theorem, we state an impor- 
tant lemma about the spectrum of Toplitz operators (explaining where the spectrum of 
the truncated matrices T n eventually concentrates). 

Lemma 5 (Spectrum of Toplitz operators and truncated Toplitz matrices | 53 |). 

Suppose A G PC (T) and let T(A) be the infinite Toplitz matrix associated to the symbol 
A. Assume U C (D is an open set and assume that the spectrum cr(T(A)) is a subset 
ofU. Then there exists an index no G N such that n > no implies a(T n (X)) C U. If 
A is continuous the spectrum of T '(A) follows entirely from geometric properties of the 
symbol, 

<j(T(A)) = range(A) U {x G € | wind(A, x) ± } . (72) 



Szego's Theorem holds for symbols in the Krein algebra l53l : 
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Theorem 2 (A second order Szego theorem (5TJ|53)). Let /i G W n B% be a symbol 
and let T n G (D nxn be the family of associated Toplitz matrices. Let Q C (D be an 
open subset that contains the spectrum ofT(fi). For an analytic function h : Q —± <D 
the following trace formula holds, 



with 



tr (h(T n )) = nGhW + ^(aO + o(l) (73) 

G h ( f i) = ^j\( f i(e w ))d6, (74) 

^(/i) = / h(X) ^- log £7(/i - A)dA, (75) 

oo 

J5(/x) = exp^fc(log/i) fe (log/i)_ fe , (76) 
fc=i 



where 



(log M ) fc :=^JJe- ikd log fi(e w )d6 (77) 
w ^ Fourier transform of \og(/i). 



Theorem 3 (Alternative expression for the second order term |52 |). Let fi G W fi 

i 

£> 2 2 Z?£ absolutely continuous and let T n G <D nxn be the family of associated Toplitz 
matrices. Let Q C (C be an open subset that contains the spectrum ofT(p) and let 
h : Q C 1 be analytic on Q. Then the second order term Eh(fi) from Theorem^can 
be written as 



-i f>7T /»7T 

= 7^ E / d( t> de 



/z(e if >) - /u(e^) 



fc=i 

x sin(fc(0-</>))( M '(</>)- M '(0)). 

Again, the derivative is to be read as in Eq. ( [37] ). 

Proof of Lemma^ We are now in the position to prove Lemma[T] We can now collect 
the results from the previous sections. By the results from Subsection |4.2| specifically 
Lemma [T] we know that for any subset of sites whose cardinality is much larger than 
/?, we can work with infinite truncated Toplitz matrices and their symbols. The asymp- 
totic expression for the entropy of the full system can be obtained directly from the 
continuum approximation. Remember that L4| = \qN] and subsystem B is of size 
\B\ = N — \A\, where q G (0, 1). Theorem [5] states that the block entropies are asymp- 
totically equal to 

S a {pA)=tr(s a (X\ A ))=qNG aa {\) + E aot (\)+o{l) (78) 

and 

SaiPB) = tr{s a {X\ B )) = (1 - q)NG aa (\)+E aa (\) + o(l). (79) 

We now turn to the computation of the entropy S a (p) of the entire system. This is 
subtle, and one cannot employ the same formula for the larger system of A U B: Lemma 
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[T] would no longer be valid, and the boundary conditions would not be respected. But 
we can still find an asymptotically exact expression. The spectrum of X is given by 



Hence, 



{(foe)(e 2 ^ N ):k = l,...,N). (80) 
S(P)= £ (s a ofoe)(e 2 ^ N ). (81) 



Using Lemma [2] again, we find that this expression can be replaced by 

S(p)=NG s (X) + o(l). (82) 
Combining the terms, the mutual information is asymptotically equal to 

I a (A:B) = 2E Sa (\) + o(l). (83) 

Since the symbol fi := foe is analytic (and therefore absolutely continuous) Theorem[3] 
can be applied resulting in the desired expression 

00 ro (A(e")) - , a (A(e*)) 



1 f u r u 

I a{ A:B) = — 2 Y, dj> 

Z7T U-A J27T J2-K 



x sin (k(0 - </>)) (A'(e i0 ) - X'(e w )) + o(l). 
which completes the proof. □ 

4.4. Simplifying Widom's second order expression. We have almost proven the main 



theorem, but we still need to eliminate the infinite sum in k in Eq. ( [38] ). The sum in k is 
rather awkward; since the integrand is more and more oscillatory for larger and larger k, 
it is a priori far from clear where one may truncate the sum in order to arrive at a reliable 
result. In this way, the expression cannot be easily computed, not even numerically, and 
practically only for low temperatures. We will hence go a technical step further and will 
prove the validity of Theorem [T] the main result for one-dimensional systems. A first 
step in this direction is the following observation. 

Lemma 6 (Function kernel). For any n e N and any 0, e [0, 2ir), 

ST • n m\ lcos(0/2) -cos((n + 1/2)0) „ f 

= 2 sin( 0/2) = : K "W' (84) 



Proof. We first note that the sum over sine functions encountered in Eq. ( 84 ) can be 
brought into a closed form reminiscent of the Dirichlet kernel. Indeed, if we truncate 
the sum at n, then it can be expressed in closed form as 



J2 sin(^) = i ]T (e ifc ^ - e" ifc *) (85) 

l<k<n 
l<fe<n 

I (\ — e O+l)i</> \ _ e -(n+l)i0 N 



2i 

Kfc<n Kfc< 



2i\ 1 - e~ 1( P 

1 cos (0/2) - cos ((n + 1/2)0) 

2 sin (0/2) " 
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where the last line follows from elementary trigonometric identities. □ 

That is to say, the mutual information can be re-expressed in terms of the kernel K n 



as 



I a (A:B) 



1 



lim 

n— )-co 



/TT rTT 
# / 
-TT J —TT 



dO- 



(\(e ie )) - s a (A(e^)) 



A(e ie ) - A(e^) 
x (A'(e^) - A'(e w )) K n ((9 - </>)) + o(l). 



(86) 



We can further simplify the expression by invoking results from the theory of highly 
oscillatory integrals. Recall the subsequent fundamental lemma. 

Lemma 7 (Riemann-Lebesgue). Let f G Z>i([a, 6], (D), let g G Loo(R, (D) be periodic 

(with period T) and assume J Q T g(x)dx = (the probably most important special 
instances of such functions g are the trigonometric functions sin and cos) then 



lim / f(t)g(nt)dt = 0. 



(87) 



Proof If / G Cc?([a, 6], (D) (hence is a compactly supported smooth function from 
[a, b] to (D), then partial integration gives 

rb -j />6 rnt 

/ f(t)g(nt)dt = — f'(t) g(s)dsdt. (88) 

for every 



Note that /'(t) is bounded by assumption and that | J Q n #(s)ds| < T 
tel. Hence, 

^6 



lim / f(t)g(nt)dt = 0. 



(89) 



For general / G Li([a, 6], (D), note that the smooth, compactly supported functions are 
dense in Li([a, 6], ©). Hence, for an arbitrary e > there exists /2 G Cc°°([&, 6], C) 
with 



11/2- /111 < 



(90) 



Therefore, 

r b 



lim 

n— >-oo 



/ f{t)g(nt)dt < lim ( / (/(t) 



h{i))g(nt)dt 



+ 



I f 2 (t)g(nt)dt 

J a 



< e 



(91) 

by validity of the statement for compactly supported, smooth functions. Since e > 
was arbitrary, the result also holds true for general / G Li([o,fc],C). □ 

This well-known lemma simplifies the understanding of the limiting behaviour of 
these kernels K n action on sufficiently regular test functions. 

Theorem 4 (Distributional convergence of the kernels K n ). Let h G Za(T, ©) £e 
locally Lipschitz continuous at e l6> G T f/ien 



lim r h{?*)K n {<j>-o)d4=\ r 

n ^°°J-TT 2 J-TT 



w ft(e^) - ft(e ig ) 
tan((0-6>)/2) < 



(92) 
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Proof. Without any loss of generality let — (otherwise take 0' := — 6). Assume 
h(l) = first. By Lipschitz continuity of h, there exists some e > and some constant 
L > such that |/i(e 1( ^)| < L\(j>\, whenever < e. By continuity, there exists some 
constant M > such that 



tan(0/2) 



< M 



(93) 



for 6 G 



7r, 7rl. Therefore 



h (e i( ^) 



tan (0/2) 



< 



X(- £ , £ )0)M • L + X[- 7r , 7r ]\(- £ , £ )(0) 



|ft(e^)| 
tan(e/2)' 



(94) 



where \a is the characteristic function of the (measurable) set A. Define h((j>) := 
h (e 1( ^) / tan (0/2). By the previous estimate h G Li([ — 7r, 7r], (D) and the Riemann- 
Lebesgue Lemma can be applied. Applying the Riemann-Lebesgue Lemma twice and 
using elementary angle-sum identities for trigonometric functions yields 

nm r 1 cos W 2)- cos ((n + 1/2)0) 1 r 1 - cos^) 

n^°°J- n x '2 sin(0/2) r 2 ]_„ v ; tan(</>/2) r 

i |" Me^)(l-cos(n^))c 



lim 

n— )-oo 2 



j — 7r tan(<) 



tan(0/2) 

If /(l) ^ write /i(e^) = (/i(e^) - /i(l)) + apply the former result to the first 
summand and note that 



(95) 



71 cos(0/2) - cos ((n + 1/2)0) 
sin(0/2) 

for every n G N by anti-symmetry of the integrand. □ 



d(j) = 0, 



(96) 



We need an extension to double integrals, a proof of which follows the same line of 
reasoning. 

Theorem 5 (Distributional convergence of the kernels K n for double integrals). Let 

h G Li(T 2 , (D) be Lipschitz continuous in a neighborhood of the diagonal 



D T 2 := {(x,y) eT 2 \x = y} 



then 



lim r ^ h(e i ^e w )K n ^-0)d^dO=^ f f 

J — 7T J — 7T J — 7T • ' —71 



h(e i( t>,e w )-h(e i0 ,e w ) 



tan 



(97) 
(98) 
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As a consequence, Eq. ( [86] ) can be rewritten in the following way: 

V 7 4ir 2 J_ n ^J_ n A(e^) - A(e^) tan ((</>- 6>)/2) W 

(99) 

Let us conclude this subsection by commenting on the second order term in the trace 
formula of Theorem [3] Recently, another related approach aiming at a simplification 
of Widom's formula has been discussed in the literature. The main result of a recent 
masters' thesis 1 56 1 is a derivation of a second order trace formula from asymptotic 
inverses and appropriate factorisations of the symbol. The presented final expression 
of the second order term does not contain k-sum as well and finally only involves the 
calculation of an appropriate double integral - a result similar to ours. However, the 
given double integral always involves complex valued integration whereas our formula 
contains only real valued integrands right from the beginning (whenever the symbol is 
real- valued and the analytic function restricted to the reals is also real- valued as usually 
satisfied in problems of practical relevance). 

Indeed, the idea to approach the trace formula starting with the inverse function 
rather then the usual logarithm (the term E(n) from Szego's theorem is just the second 
order term of the logarithm) has many didactical and computational advantages in our 
opinion. First of all, this computationally inappropriate and practically incomputable 
k— sum in Eh(fj) is missing right from the beginning - note that for larger and larger k, 
the integrand becomes more and more oscillatory in a fashion that is hard to grasp. Fur- 
thermore, there exist well-known asymptotic expansions for asymptotic inverses even 
for higher order terms (which are very useful for approximate solutions of systems of 
linear equation with coefficient matrix being a Toplitz matrix for example). Last but 
not least the authors of Ref . l56l were able to derive second order formulas for block 
Toplitz matrices (which require - unlike the final result in the one-dimensional setting 
- an explicit factorisation of Al — /i where \i is the matrix- valued symbol and A is an 
arbitrary complex number on the integration contour). 



5. Temperature dependence 



5.7. General remarks. Our main result, Theorem 4.1 gives an asymptotically exact 
expression for the mutual information of two neighbouring blocks of fermions in one 
dimension. It has been shown that for a given inverse temperature f3 > 0, 

I a (A : B) = I^A : B) + o(l) (100) 

in the system size N. The result given constitutes an easily computable expression that 
paves the way for studying a number of physically meaningful regimes. At this point, 
we would like to pause for a moment, however, and would like to come back to one of 
the questions posed in the introduction, namely of the possible asymptotic behaviour 
of I a ,oo(A : B) for large and small inverse temperatures. As we have seen above, the 
general bound of Ref. l42l following from the extremality of the free energy suggests 
that the mutual information should scale like fj for large temperatures, 

I ajOO (A:B)=0(p). (101) 

One might wonder whether this bound actually gives the proper asymptotic scaling on 
the temperature. Conformal field theory actually suggests a behaviour which is logarith- 
mic rather than linearly in the inverse temperature EJ [27]|. In this section, we corrob- 
orate the prediction from conformal field theory by showing that the low temperature 
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Critical phase (with b = 1, a = 1) Non-critical phase (with b = 1, a = 4) 




High temperature asymptotics 
(where fli(a,6) := c*6 2 /(2 log(2))) 

J li00 (A : 5)//3 2 /i,oc(A : B)/l? 




0.2- 0.2- 



0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 0.5 

Low temperature asymptotics 
(where R 2 (a) := a&/(ir* log(2)) min { (min{ 2 J 1})2 , ^ } , ife(a, b) := - (|a| - 26)) 

/ lj00 (A:S)/ln(/3) ln(/ li00 (A;B))/^ 




Fig. 2. Temperature dependence of the von-Neumann mutual information for various parameters of the 
fermionic instance of the XX model. 
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asymptotics are given by 

I a , 00 (A:B) = 0(log(f3)). (102) 

This has a quite remarkable consequence: In order to see features of the criticality of 
the ground state, one has to go to extremely low temperatures. This dependence is also 
convincingly depicted in Fig.[2j where a logarithmic scale has been chosen - otherwise, 
signatures of ground state features could hardly be detected. Only at extraordinarily low 
temperatures, the familiar logarithmic divergence in the system size of the sub-block 
chosen becomes visible for a reasonably sized subsystem. 



5.2. Analysis of the XX model Needless to say, the low temperature asymptotics de- 
pends on the choice of the model parameters a and b. Specifically, the scaling of the 
mutual information reflects signatures of the zero temperature quantum phase transi- 
tion of this model taking place at \a\ = 2b. Whenever \a\ < 26 the model is critical in 
its ground state, i.e., there is no energy gap. It is known that the entanglement entropy 
then exhibits a logarithmic divergence - signatures of that are also seen in the mutual 
information at small but non-zero temperature. If the model is gapped, i.e., whenever 
\a\ > 2b, the ground state is the vacuum or the fully occupied state, depending on the 
sign of a. As expected, the mutual information converges exponentially quickly to zero 
in the non-critical phase at a rate proportional to the energy gap. The main insights for 
this model are summarised in the following theorem: 

Theorem 6 (Temperature dependence). The (a—) mutual information of the thermal 
state of the fermionic XX model with parameters a G R and b > satisfies 

Inn f °."% ; *> = (103) 
(3^0 p 2 21og(2) 

Whenever \a\ < 2b, then 

r la, oo\ A : B) rif\A\ 

limsup — — (104) 

^oo log P 



2a 

< o mm 

7T 2 log(2) (min{l,a}) 



( 27 4 \ 

2 111111 \2(min{l,a})' ]l^a~\ J * 



Whenever \a\ > 2b, then for any n < (\a\ — 2b) • min{a, 1} 



lim exp(K,p)I a oo(A : B) = 0. (105) 



The proof of this statement is rather involved and requires a number of techniques 
developed in lemmas: Hence, for better readability of the main text, it will be presented 
in the appendix in Subsection [84 



Area laws for thermal free fermions 



23 



6. Free fermonic models on the torus 

The results established above readily apply to the situation of a bi-sected, higher-dimen- 
sional fermionic lattice system on the torus. A quite similar strategy has already been 
exploited, e.g., in Ref. l25ll . Using appropriate discrete Fourier transforms, one can 
disentangle the constituents from each other with respect to all but one dimensions. The 
exception is the dimension in which the two regions labeled A and B are singled out, 
see Fig. [3] In this way, one arrives at a collection of suitably modulated and altered 
one-dimensional problems, to which the above statements apply. In this section, we 
highlight the results obtained in this manner. An important application of this is the 
computation of the mutual information in higher- dimensional tight binding models. 



6.1. Geometry of the problem. Let us for simplicity consider in D dimensions the ge- 
ometry of slabs L = (7Lm) D ~ x x ^at, for suitable N and M. So along one dimension, 
we have as before TV sites, whereas the system size with respect to the other dimensions 
is M, see Fig. [3] The index set of all sites can be taken to be 

I = 1' x J, (106) 
/ / = {0,...,M-l} xD - 1 , (107) 
J = {0, . . . , AT — 1}. (108) 

Vectors of indices i G / C 7L D will be regarded as modulo N and M, respectively. Let 
w : 7L D — > L be defined as 

w(i) = w{i\, . . . , %d) = (iimodM, . . . , i £>_imodAf , z^modA/'). (109) 

This function simply projects arbitrary indices from 7L D onto the lattice L. We again al- 
low for arbitrary finite-ranged interactions (the generalization to exponentially decaying 
interactions is straightforward but omitted). That is to say, for z, j G /, the Hamiltonian 
tensor takes the form 

Vj=d w{l _ 3) . (110) 
Then the Hamiltonian can be written in the following way 

H -= E vtffi = £ v^fl c) f {b4) . (in) 

ijel a,beI';c,deJ 

We have singled our the special spatial dimension for which we consider the bi-partite 
cut. We take the parts A and B to be 

A = / / x{0,l,...,^A-l-l}, (112) 
B = V x {\qNl\qN] + 1, . . . , N - 1} , (113) 

with q as before. We will also see that the previously found results still apply. 
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Fig. 3. The geometry of the free fermionic model on the torus. Along a singled out dimension, the system 
consists of TV sites, while the other dimensions with the topology of a torus embody M sites. The distin- 
guished region is again referred to as A, its complement is B. 



6.2. Scaling of the mutual information on the torus. We can now consider Fourier trans- 
forms with respect to the index set I\ while leaving the index set J invariant. For the 
discrete Fourier transform, we choose 

E/ M := 7 ±= exp (J^k • l) for k, I e I'. (1 14) 

It can then be shown that 

{ M v ' v 3 >i,i v > ) (k,k ),(iM 

= (5^exp(^t-z)d (t>feo _, o) J 5 kJ . (115) 



\tei' 



=4° , 

k ,l 



That is to say, the reduced entropies of A and B as well as the global entropy can be 
computed as if one had uncoupled one-dimensional systems, only that the coupling is 



modulated in a way described by the new couplings in Eq. ( |115| ). Denote for k G I' the 
resulting expression of the symbol computed from 

Afe = /o£: fe , (116) 



where Sk is the expression as defined in Eq. (20), with d replaced by d^ for each 
k. In this way, we also arrive at the asymptotically exact expression for the mutual 
information in the limit of large N, recovering the logarithmic divergence of the mutual 
information. 



Theorem 7 (Mutual information on the torus). For any inverse temperature > 0, 
any M G N, and for any a G [0, oo), the mutual information is given by the, in N 
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asymptotically exact, expression 

^ :B ) = iE [ d( t> f de 

Sa(A fc (e w )) - a a (\ k (eH>)) \' k {4>) - X' k (0) 



x 



\ k (e«>) - \ k (e^) tan ((0-^/2) 



o(l). (117) 



In particular, in this way, we find for tight binding models on the torus in any dimension 
D that 

I(A : B) = 0{M D ~ 1 log(/?)) (118) 

for large inverse temperatures /?: The logarithmic divergence in the inverse temperature 
remains, while an additional term relating to an area law in D dimensions emerges. 



7. Outlook 

7.1. Conformal field theory and entanglement spectra. In this section, we compare our 
results with predictions from conformal field theory. Indeed, our findings can be seen as 
a confirmation of the predictions resulting from the conformal transformation relating 
finite systems at zero temperature to infinite systems at finite temperature. In our result, 
we provide system specific qualifiers of the asymptotic limits, and our main theorem 
constitutes a fully rigorous result applicable to a large class of models. Also, higher- 
dimensional systems on tori can be captured with the methods presented here. Still, it is 
interesting to see that the behaviour of the mutual information, scaling asymptotically 
as the logarithm in the inverse temperature, can also immediately be suggested from 
the analysis of conformal field theories in 1 + 1 dimensions. This is an observation 
similar to the one made for ground state properties of the XX model in the ground state 
(TOlES): Here, the rigorous expression given relates to and confirms the formula 0161 

SU = £log- + ci, (119) 
3 a 

in a connection already conjectured in Ref. |8|. Here, c\ > is a constant and a is 
the lattice spacing, while c is the conformal charge. The results on entanglement en- 
tropies in the XX models were proven using Fisher-Hartwig type methods for Toplitz 
determinants lfT0ll26l . 

The main idea in the conformal analysis is that there exists a connection between 
correlation properties of finite systems at zero temperature and infinite systems at fi- 
nite temperature. Therefore, assume that at zero temperature, the entanglement entropy 
scales as 

Sa = fil) (120) 

for some function / : R + — >> R + . Then, at large inverse temperatures fj > 0, we 
consider the conformal mapping z \-> e 27TZ / for v > 0, to arrive at the expression 

G3 

S A = f( V J>^(*)). (121) 



7T \ V 

Combined with Eq. (|120|), one hence gets for > I2ll27l 



S A = C -\og(^smh( 2 -^) ) + C3 , (122) 



6 I ™ I 
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for c > 3 > 0. The above statement becomes meaningful in the limit when I ^> f3. 
Asymptotically, Sa is well approximated by ttcI/(3/3) + cs [2], but for our purposes, 
we also need a first order correction, namely 

s -H( los GQ- log(2,+ T) +C3 - < 12 » 

In this way, one recovers the expression for the mutual information 

I(A:B)*^log(^-), (124) 
3 \7ra J 

again exhibiting a logarithmic dependence in /?, as (j —> oo. The analysis of the finite 
temperature case axing arguments from conformal field theory has recently been carried 
out in Refs. (58] [59). ^ s men ti° ne cl earlier, the criticality of the ground state is not felt 
by the mutual information, unless one is at extremely low temperatures. In turn, the 
logarithmic scaling in the size of a subregion is turned into a logarithmic dependence 
on the inverse temperature. 

It is important to note that because the formalism developed here allows for the 
computation of all Renyi entropies, the tools developed in this paper also allow for 
the study spectra of reductions; i.e. entanglement spectra lt36l44TTl . compare also Refs. 
HJ |3 . Such entanglement spectra have turned out to provide a powerful tool when 
characterizing topological order, relating bulk to boundary theories, and discussing the 
possibility of approximating states with suitable tensor network states such as matrix 
product states and operators and higher-dimensional analogues. 



7.2. Remarks on open free fermionic many-body systems. Let us finally briefly mention 
that the formalism developed here in principle also allows for the study of open quantum 
systems, undergoing dissipation and quantum noise. An extensive discussion of this 
topic is beyond the scope of the present article. In the mindset developed here, one can 
consider Liouvillians capturing fermionic open quantum systems of the form 

C(p) = i[H, p]+Yl i L 3P L ) ~ \{L]L^p}) (125) 

with, as before, 

and the Lj are operators linear in the fermionic operators, supported on finitely many 
sites only, which are all the same, except that they are all translates of each other. The 
Lindblad operators Lj hence act locally in the same way as the Hamiltonian terms act 
locally. The equations of motions then become 

j t P{t) = C{p). (127) 

A state p is called a steady state if C(p) = 0. For fully translationally invariant systems, 
one can again define a symbol OHS), and the methods developed here are applicable. 
Details of such an approach will be pursued elsewhere. Such an idea seems particularly 
timely in the light of the observation that open system dynamics should not only be 
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viewed as a source of quantum noise added to the system. But that open system dy- 
namics and dissipation can also be beneficial: Indeed, dissipative dynamics allows for 
dissipative instances of quantum information processing [46 -48 ], exhibiting an inter- 
esting way of protection. Also, steady states can readily exhibit a number of exciting 
properties: they can be entangled, show phenomena of noise-driven criticality, or even 
exhibit topological order I44H461 . The latter is specifically true for free fermionic mod- 
els, an arena for which the machinery developed here should be fruitful. 

7.3. Summary. In this work, we have introduced a formalism allowing us to prove 
the validity of rigorous expressions for the Renyi mutual information of Gibbs states 
of translationally invariant quasi-free fermionic models. The expressions obtained are 
asymptotically exact; the bounds given are also exponentially tighter than those derived 
from the extremality of the Gibbs state with respect to the free energy. Special em- 
phasis has been put onto the technical development of novel methods of dealing with 
Toplitz matrices, in particular a new and useful instance of a second order expression 
for smooth symbols. These tools, as well as the approximation results introduced here, 
are expected to be widely applicable and highly useful in various contexts. It is the hope 
that this work inspires further entanglement-related studies of mixed fermionic quantum 
states, arising in the context of describing both closed and open quantum many -body 
systems. 
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8. Appendix 

8.1. Simple and useful upper and lower bounds to the mutual information. In this 
subsection, we will derive upper and lower bounds to the quantum mutual informa- 
tion, evaluated for Gaussian fermionic states. These bounds are generally useful when 
bounding mutual information expressions. Moreover, they are additive, so share this im- 
portant feature with the actual mutual information. We expect these bounds to be useful 
also in other contexts, different from the study pursued here. What is more, it should be 
clear that one can immediately also formulate a bosonic variant of the bounds presented 
here, using the same strategy of proof. We consider covariance matrices 7 G R 2nx2n 
of the form 

i=(-xo)> < 128 > 



with X = X T . We write 



(129) 



and denote with P the pinching l55l 



X A 

x B 



X A ^X B (130) 



of X. This is the covariance matrix of the tensor product of the two reduced states, as a 
moment of thought reveals. In these terms, the quantum mutual information can hence 
be written as 

I(A : B) = tr(s(P)) - tr(s(X)), (131) 

so as a difference of two trace functions. Let us now consider a quadratic function 
I : [-1,1] [0,1] defined as 

l(x) = (l-x 2 )/2. (132) 

We will see that the corresponding difference of trace functions is a convenient, com- 
putable lower bound of the mutual information with reasonable properties. We hence 
relate the mutual information to purities only (or, for that matter, to 2-Renyi mutual 
informations J] 



2 It may be important to state that since these quantities only referring to local and global purities are 
easier to measure, these bounds are expected to be useful even in an experimental context. 



30 



H. Bernigau, M. J. Kastoryano, and J. Eisert 



Lemma 8 (Lower bound to the mutual information). For any X = X\ 

tr(s(P)) - tr(s(X)) > tr(/(P)) - tr(/(X)). (133) 
Moreover, the bound is additive, in that 

tr(Z(P)) = tv(l(X A )) + tr(l(X B )). (134) 
Proof. To start with, the additivity immediately follows from the definition. We have 

tr(/(P)) = itr(l 2n -(X A ©X B ) 2 ) 

= \tr(t 2n -(X\®Xl)) 

= tr(Z(X A ))+tr(Z(X B )). (135) 

This lower bound as such can be shown using a result of the theory of convex trace 
functions called Peierls' inequality l55l . The function s — I is a concave function (it is 
not required that it is operator concave). What is more, we can without loss of generality 
assume that Xa and Xb are both diagonal. If they are not, there exist real orthogonal 
matrices Oa,Ob such that (Oa © Ob)P(Oa © Ob) t is diagonal, without altering the 
value of the mutual information or the bound. Using the concavity of g := s — I, we 
find 

J29(X jd )>tr(g(X)). (136) 

3 

But since Xa and Xb are diagonal, we have 

tr(g(P)) > tr(g(X)) = tr(s(X)) - tr(l{X)) (137) 
which is what we intended to show. □ 

Note the similarity of this bound for the mutual information with the bound presented 
in Ref. [14] for the von-Neumann entropy. It is quite remarkable that this bound can 
even be used for differences of entropies, not only of entropies as such. 

We now turn to upper bounds. Let us define the function u : [—1,1] — » [0, 1] as 

tt(l) = i) (1 " l2)1/2 (138) 

Again, one can prove that the upper bound derived from this matrix function, tr(u(P)) — 
tr(u(X)) shares the above additivity property. Again, one can prove in the same way 
as above that one encounters an upper bound. Of course, the function u — s is no longer 
concave, but convex, so I — s is concave. This bound is not a quadratic bound any more. 
Using again Peierls' inequality, one arrives at the following result. 

Lemma 9 (Upper bound to the mutual information). For any X = X\ 

tr(s(P)) - trO(X)) < tr(u(P)) - tr(u(X)). (139) 
The bound is additive, 

tr(/(P)) = tr(u(X A )) + tr(u(X B )). (140) 
This upper bound on the mutual information is a very useful bound in its own right. 
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8.2. Derivatives of entropy functions. For convenience, we present the derivatives of 
the entropy functions used here. 

Lemma 10 (Derivatives of the entropy functions). The derivatives of the entropy 
functions are given by 

Sa(X > (l-a)log(2) (l +x f + (l- x f ' U4U 



a ( (1 + x) a ~ 2 + (1 - xf' 2 



log(2) ^ (1 + x) a + (1 - x) a 
a f (l + x)"- 1 -(l-x)^ 1 " ^ 



l-a\ (1 + x) a + (1 - x) a 



(142) 



= ilog 2 (^), (143) 

^-^(ih + ih)- (144) 



8.3. Second order theorem for convex functions ofToplitz operators. In this subsection, 
we elaborate on second order theorems for convex functions of self-adjoint Toplitz 
operators. It builds upon the analysis of the main text, but is not directly required in any 
of the above proofs. For the main statement of this subsection, we will need Peierls' 
inequality in the following form. 

Lemma 11 (Peierl's inequality). Let A = be a Hermitian N x N matrix and let 
U be a unitary N x N matrix. Then, for any function f that is convex on some open 
interval I with a A G /, one finds that 

J2 f({UAU^)..)<tvf(A). (145) 

l<i<AT 

The following statement holds. 

i 

Theorem 8 (Self-adjoint Toplitz operators). Let \i G W ^\B^ be a symbol. Assume 
further that range(/i) C Rfor equivalently assume that T(/i) = T(p)^). Then the 
associated family ofToplitz matrices, T n £ (C nxn is Hermitian. Let J CRk an open 
interval with range (/i) C J and let h : J — » R be real- analytic. Then there exists 
an analytic extension h of h to a small complex neighborhood of J of range (a) with 
h |range(» = ^ |range(/x) (compare Lemma\2v and the trace formula ofTheorem\2\holds 
true for this extension. Moreover ifh is convex, then 

E h (jj) < (146) 

Proof. The only non-obvious statement is 



E h (tJL) < 



(147) 
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Fix e > and no G N large enough such that 

trft(T n ) = nG h (ii)+E h (fi)+r n (148) 

with \r n \ < e/3 for every n > uq. Since T n is Hermitian we have r n G R. Choose 
unitary matrices C/ that diagonalize T no . Then 

2h(T no )= h((U®UT 2no Ul®U^..)=2n G h (ri + 2E h ( f i) + 2r no . 

l<i<2n 

(149) 

Hence, by Peierl's inequality: 

2nG h (n) + 2E h ( t i) + 2r no < tr (ft (T 2no )) = 2r^G h {ii) + ^(/i) + r 2no . (150) 
Therefore: 

#/i(aO < r 2no -2r no < e. (151) 
Since e > was arbitrary this implies 

E h (fi) < 0. (152) 

□ 

Since obviously E hl+h2 (^) = E hl (fi) + ^ 2 (/i) and E X h(v) = ^Eh(v) for A G © 
this implies that E^ x (/x) < i^ 2 (/i) whenever fti — ft 2 is convex. If both fti and ft 2 are 
differentiable twice this holds true if and only if ft" < ft 2 '- 

8.4. Proof of Theorem^ 

Proof. We start with proving the part of the theorem relating to the high temperature 



dependence, Eq. ( |103| )). For that purpose, we formulate the explicit asymptotic expres- 
sion for the mutual information of a thermal state of the fermionic instance of the XX 
model. Define 

e^):=/8(|+6cos(0)). (153) 

Then 

4(0) := -06sin(0). (154) 

The symbol takes the form A(e 1 *) = — tanh(e J g(<^)). In these terms, we can express the 
mutual information as 

Ia t00 (A : B) = ^ J* J d0d4>fp(8,4>)tp(0,<P), (155) 

where we have decomposed the integrand into the functions 

r ,n A \. a a (tanh(e^(6>))) - s a (tanh(e^(0))) 

■ m 0) - toipjH^ (156) 
\ 

cosh 2 cosh 2 ( e/3 (0))/ ' 

^(0,0) -tan" 1 (157) 
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The high temperature (small (3) limit can be calculated simply by expanding the func- 
tions in the integrand, to leading non-zero order in /?, which yields 

MO, <t>) = (a + b (cos 9 + cos <t>)) (sin 9 - sin </>) + 0(/3 3 ). (158) 

Now, noting that 

/ d(9d0 (a + 6 (cos (9 + cos 0)) (sin (9 - sin <j>) tan -1 (((9 - <j>)/2) = 46tt 2 , 

-7T J — 7T 

(159) 

we find that 

lim I a MA :B) = a* 
(3^o f3 2 21og(2) 

□ 

The second part of the theorem, relating to the low temperature asymptotics (i.e., 
the asymptotics for large /?), turns out to be significantly more involved. In a first step, 
we will cast the integral into a more appropriate form. By an elementary argument 
involving the symmetry of the problem, one finds that 

I atOQ (A:B) = ^(I 1 +I 2 ), (161) 

where 



and 



h = r r dBd4Mw)Hw)> ( i62 > 

Jo Jo 

nTT pO 

h= / d0dMp(0^)tp(6^), (163) 

JO J -TV 

In this expression, we have made use of the same notation as in the above proof of the 
low temperature behavior. Note that both cos functions are monotonous functions of 
there argument in both regions of integration. The substitution 

x := cos(0), (164) 
y := cos(0). (165) 

and an appropriate rearrangement of factors then gives 

h=fib J J dxdyfp ill (x,y)(fp tl 2(x,y) + fp t i3(x,y)), (166) 

where 

, / x s a (taiih(ep(x))) - ^(tanh^ (?/))) 

jp t ii(x,y) := , (167) 

x-y 

f , = sech 2 (e>(x)) (1 - x 2 ) - sech 2 ^)) (l - y 2 ) 1 

W 2 ^'^' ^ tanh(^(x))-tanh(e /3 ( 2 /)) ^ (1 - x 2 ) 1 /2(i _ ^2)1/2 ' 

(168) 

f „\ — sech 2 (e^(x)) -sech 2 (^fa)) 

tanh(e /3 (x)) - tanh(e /3 (^/)) 
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In this equation, we have defined 

ep(x) := P{a/2 + bx). (170) 
An analogous computation for I 2 yields: 

h = /3b J J dxdyf^n(x,y) (fp,i2(x, y) ~ fp, is(x,y)) . (171) 
We need the following elementary but powerful lemma. 

Lemma 12 (Bounds on difference quotients). Let f,g,h,kbe differentiable functions 
on some open interval ICR. Assume 

- s'(z) > for any s G {/, g, h, k} and any z E I, 

- sup {s'(z) \z G /} < oo for any s G {/, g, h, k}, 

- f'{z) < h'{z) andk'(z) < g' \z) for all z G L 

Then, for any pair of real numbers x,y <E I with x ^ y, 

< f(x)-f(y) < h(x)-h(y) ^ 
- g(x) -g(y) ~ k(x) - k(y)' 
Proof By the mean value theorem there exists some z £ I with x < z < y such that 

(h-f)( X )-(h-f)(y) = m _ fl{z) > Q (1?3) 
x y 

Hence 

m^m > m-m > 

x — y x — y 

where the last inequality follows from f > 0. By the same argument 

g(%) - g(y) > Hx) - kjy) > Q 

x — y ~ x — y 
multiplying both inequalities and rearranging terms yields the desired lemma. □ 
Moreover, we will use the following integral. 

Lemma 13 (Integral over exponential difference quotient). For any a\ < a 2 and 

h < b 2 , 

Aa u a 2 MM : = / / dxdy (176) 

J ai Jb! v~ x 

= F(b u a u a 2 ) + F(a u b u b 2 ) - F(b 2 ,a u a 2 ) - F(a 2 , b u b 2 ). 



Here 
with 

where 



F(a, b, c) := exp(— a) F(b — a) — F(c — a) 



(177) 



F(x) := Ei[-x] - log \x\ =\e + J2 (178) 

k>l 

Ei(x) := f ^ (179) 
is the exponential integral function and denotes the Euler Gamma constant. 
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Q 5 

Qe,5 



-0.5 ■ 



Qr,5 



-1 -0.5 0.5 1 
Fig. 4. The different regions of integration for a = b and S = 0.5. 

Having established these two lemmas we proceed with the proof of Theorem [6j 
Proof. (Of Theorem |6| Eq. ( |104| ) and Eq. ( |105| )). Consider the open interval 

J:=(f- fc 'i + & )- (180) 



26 



2b 



< 6m 



}• 



If e I (i.e., \a\ < 2b), set m := 1 - |a|/(26). If £ J, set m := |a|/(26) - 1. For 
1 > (5 > we will split the region of integration into several sets, namely the following. 

- The bulk contribution: 

Q s :={(*,?/) E [-1,1] x [-1,1] 

If \a\ > 2b, then Q 5 := for every < (5 < 1. 

- The edge contribution: 

Q e ^:={(x,2/)G[-l,l]x[-l,l] 

If \a\ > 2b, then Q 5 = [-1, 1] x [-1, 1] 

- The remaining contribution: 



2b 



2b 



> Sm 



(181) 



(182) 



Q r ,s := [-l,l]x[-l,l]\(Q*UQ M ). 



(183) 



The various regions of integration are depicted in Fig. [4] We will prove the theorem in 
several steps: 

Step 1: The Q e j contribution of Ii and I2 decays exponentially fast. Formally, we 
will show that for any < 5 < 1 and every k < 28bm min {a, 1} 



lim exp(ft/3) / / dxdyfp ill (x,y)(fp il2 (x,y)±fp il3 (x,y)) 



0, (184) 



where we used the same notation as in Eq. ( |166| ) (note that this proves Eq. ( |105| )). 
First of all observe that 



(//3,12 ± /^,13) = |//3,ll| • |/0,12 ± . 



,13| 



(185) 



36 



H. Bernigau, M. J. Kastoryano, and J. Eisert 



Consider the factor \fp,n\ first. If y > —a/26 + 5 or x, y < —a/26 — 5 (this is 
actually true for all x, y whenever \a\ > 2b) the mean value theorem yields 



|i/3,ii | = /?& • sech 2 (ep(z)) ■ s' a (tanh (ep(z))) 
< (/3b) 2 d exp(-2 min {a, 1} /3&£ra). 



(186) 



Here z lies between x and y and Ci is some constant (any C\ > 4a/log(2) works). 
Here and in the future, we will use the symbol < to indicate that < holds for suffi- 
ciently large f3. 

If \a\ < 2b then x + a /(2b) < —5, y + a /(2b) > S has to be estimated separately. In 
this case consider the estimate 



< 



2 sup {s a (z) \ \z\ > tanh(£/?bm) } 
25 



< C 2Aa (3bexp(-2(3m5ba), (187) 



for some (possibly S— and a— dependent) constant C 2y §^ a (the factor j3b is actually 
only needed if a = 1). So asymptotically, we have that 



|//3,n| < (C 1 +C 2Aa )(^) 2 exp(-2min{a,l}^m). 
Consequently, setting C 2 '= C\ + C^,^, we obtain 

/ dxdyfp t n(x, y)(fp,i 2 (x, y) + fp,iz{x, y)) 
< C 2 (b(3) 2 exp(-2/36(Smmin{a,l})J, 



(188) 
(189) 

(190) 

(191) 
(192) 

J<(2 + e)7r 2 (193) 

for any e > 0. This together with Eq. ( |189| ) yields the identity Eq. ( |184| ) and finishes 
the first step. 

Step 2: 77ze Q r j contribution of I\ and I 2 is bounded. Formally we will show that 



where 

J= dxdy\fp :12 (x,y)\-\-\fp A3 (x,y)\. 
</[-!, l]x[-l,l] 

The generalized mean value theorem and some calculus yields 

sech 2 (ep (x)) - sech 2 (ep(y)) 



< 2 



and 



Hence 



tanh(e#(x)) — tanh(e^(x)) 
sech 2 (e>(x))(l - x 2 ) - sech 2 (e>(y))(l - y 2 ) 



tanh(e#(x)) — tanh(e^(x)) 



< 2 



2 

bp' 



lim sup /3b 



dxdyfp tll (x, y) (fp,i 2 (x, y) ± fp t iz{x, y)) < 00. (194) 



Without any loss of generality we can assume that 5 < 1/4. We will will split Q r j 
into the sets 



•,5,1 



Q r ,5n{(x,y) g [-1,1] x [-1,1] ||a; < 5} 



(195) 
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and 

Qr,S,2 '•— Qr,S,l \ Qr,S,l- 

We will show separately that the two integrals 



(196) 



Ji~Pb f I ^|/ W i(^^)I(|/ W 2(^2/)| + |/ai3(^2/)I) 5 (197) 



and 



J 2 := Pb 



dxdy\fp tll (x,y)\ (\fp,i2(x,y)\ + \ fp,i${x,y)\) (198) 



are bounded. In order to bound J\ , note that by the mean value theorem 



(1/^,12^, + I//3,13(Z, V)\) < 

Hence, 



1 



J 2 < Ctfb 
An appropriate linear substitution gives 

/3b / / dxdy\fn(x,y)\ 



f3bj (l-4J 2 ) 1 /2(i_ 52)1/2 
dxd2/|/n(a;,2/)| . 



2(56/3 /.<5/36 



JSbp JO 



dxdy 



s a (tanh (x)) — s a (tanh (?/)) 



x-y 



(199) 
(200) 



(201) 



(202) 



Note that the function f a : [0, oo) —> R defined as / a (:r) := s ot (tanh(x)) is mono- 
toneous decreasing. Moreover, 



2a 



(203) 



|l-a|log(2) 
for a ^ 1 and 



exp (—2 min {1, a} x) 
-fi(x) 



1 — exp(— 2 \a — l\x) 



1 



1 + exp(— 2ax) 1 + exp(— 2x) 
xexp (—2a;) 



log(2) (l + exp(-2x)) 2 ' 
Therefore, there exist constants C a ,K a > such that 

- fa( x ) < C a exp(-K a x). 

That is to say, Lemma [T2| gives 

s a (tanh (x)) - s a (tanh (j/)) C a exp(-^ a x) - exp(-« a j/) 



< 

y-x 

giving rise to the estimate 

r 25b/3 r 5(3b 



< 



y-x 



dxdy 



< 



5b(3 JO 
C a A 



s a (tanh (x)) — s a (tanh (y)) 



x - y 



(204) 



(205) 



(206) 



(207) 



2^k q ,25b(3K a ,0,5b/3K a > 
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where ^ai,a 2 ,6i,6 2 has been defined in Lemma 13 By some calculus again 

lim Asbp Ka ,25bpK a ,o,5bp Ka = log (2) • (208) 

P^oo 

Therefore, limsup^^ J\ < oo. In order to bound J2, the singularities of fis on 
the lines x = 1 and y = 1 have to be controlled. We use the bound 



\h^,v) + A,sl < (4 + ^) (^^v^.^ 

again and bound in the following way: 



(209) 



(1-£ 2 )V2(1_ 1/2)1/2 



1^ x ( x y )| < ^(tanh(e /3 (x))) + g a (tanh(e j9 (i/)) ^ 1Q) 

A short calculation shows that 

4 + — r 1 r 1 

^2<8 — /^(pb) 2 C a / zexp(-2/?aAaO 
Jo Jo 

(211) 

where C a is some constant with C a > 4a /log (2). Note that 

rl rl 11 7T 

lim sup f3 2 / / dxdyxexp(-2/3bax) T ^-77^73 ^tt^ < 0/0 m - 

/3^oo 7o 7o (1 -£ 2 ) 1/2 (1 -Z/ 2 ) 1 ' 2 2(2a6) 2 

(212) 

This follows easily from explicit integration of the y— integral, giving 

d ^T\ ^2U72 = \ (213) 

(l-y 2 ) 1 / 2 2 

and an appropriate estimate for the x— integral, namely for fixed e > note that 

1 1 

dxxex.p(— 2f3bax)- , (214) 

(1 — x 2 ) 1 ' 2 

< (!_g2)i/2 J dxxe^(-2pbax)+exp(-2pbae) J _ 1 /2 

^(2a6/3)2(L C 2)i/2+^(-^)f- 

Finally, multiply this expression by 1 , take the lim sup for fj to infinity and let e 
approach zero afterwards (bounding 1/(1 — x 2 ) 1 / 2 by 1 from below actually shows 
that Eq. \2\2) holds true with equality and that the limit superior is actually a limit). 
In any case J 2 is bounded, proving the second step. 
Step 3: The integrand fp, 11 bounds the asymptotic growth ofI\ and I2. Formally we 
will show that 

h ^ v . 4pbf Q dxdy\fp A1 (x,y)\ 
hm sup 1 — — < hmmf — — (215) 

^oo log(/3) P^oo log(P) 



Area laws for thermal free fermions 



39 



and 

h ^ r • 4fib J dxdy 1/(3^x^)1 
hm sup - — — < hm mf ^— — — (216) 

^oo log(/3) P^oo log(/3) 

for any < 5 < 1. The result for I2 is trivial when the result for I\ is settled, since 
I2 < Ii . Fix e > and note that by Step 1 and Step 2, we find that 

lim sup 7l = lim sup Il,e (217) 

^oo log(/3) /3^oo lOg(p) 

where 

h,e--=bf3 J dxdyfn(x,y)(fi2(x,y) + f 13 (x,y)). (218) 
However, for x, 2/ G Q e the mean value theorem yields the estimate 

|/i2(ar,y) + /i 3 (a:,tf)| < 2 + 2(l + (219) 

again (compare Step 1 and Step 2). This is especially true for the choice < e < S, 
and therefore 



lim sup - — ^r— < liminf (2 + 2- — 



fibf Qe dxdy\fp A1 (x,y)\ 
log(/3) " T-V~ V~ ' "l - log(/3) 

/ 1 \ Pbf dxdy\fp in (x,y)\ 

< liminf 2 + 2- JQ < . (220) 

/3^oo v i-e 2 y log(^) 

In the last expression e can be sent to zero, proving the upper bound, Eq. ( |215| ). 
Step 4: Derive a logarithmic bound for I\ and show that I2 increases sub -logarithmically. 
We will show that 

h . f 108a 16a 1 

lim sup - — — — < mm < 



^oo log(/3) [log(2)min{l,a} 3 ' log(2) |1 -a|min{l,a} 2 j 

(221) 

and 

lim sup - — 2 — = 0. (222) 
^^oo log(/3) 



We use Eq. (215 ) from Step 3, with which we make a linear substitution and exploit 
the symmetry of the integrand to get 

hm sup - — — < hm sup — — , (223) 

/3^oo log(p) p^oo 10g(/3) 



where 



and 



j, 1= r rw°w-^> ,224, 



^rrw-" 1 ^-" 1 <22 5 ) 

Jo Jo x + y 
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Here, f a are the functions introduced in Step 2. Fix < e < min {a, 1}. By Eq. 
( |205| ), one finds 

— fa'(x) < C a exp(— K a x). (226) 
A look at Eq. ( |205| ) and Eq. ( |204| ) in turn reveals that 

2a 

: = i — 77^? := 2(min {a, 1} — e) 



log(2)e 

is a possible choice. Proceeding as in Step 3 yields 

i • Jll . i . C a Aq SbBK a ,0,5bpK Q 

lim sup - — —— < lim sup 



'^oo log(/3) ^^oo (^ a ) 2 log(/3) 
and optimizing for e gives 

r Jn . 27a 

lim SUp : r^rr" < 



log(/?) 81og(2) (min{l,a}) 
Another possibility (working for a ^ 1 only) is to use the bound 

2a 



fa ' {x) -\l-a\log(2) 
giving the asymptotic bound 



exp(— 2 min {a, 1} x), 



r Jll . OL 

lim sup - — — - < 7y . 

log(£) |l-a|log(2) (min {a, 1}) 

To bound J12, note that x + y > \x — y\ for any #,2/ > 0. Therefore, 

J12 < Jn- 



(227) 



(228) 



(229) 



(230) 



(231) 



(232) 



To show the validity of Eq. ( |222| ), we use the mean value theorem to yield the fol- 
lowing bounds 



sech 2 (e /3 (x)) - sech z (e^(y)) 



tarih(ep(x)) - tanh(e^(x)) y^/l - x 2 y/l - y 2 
for diWx^y G Q5. By the same argument, 



< 2 



1 - J 2 



x' 2 sech 2 (e(3(x)) - y 2 sedi z (ep(y)) 



2/~ 



tanh(e /3 (x)) - tanh^fV)) ^Vl - - y 2 ^ 

and therefore 



< 



26 , 2 A2 
gb ^ Z0 

l-£ 2 



I/12-/13I <2 



^ 2 + ^ ' 
l-£ 2 



(233) 



(234) 



(235) 



repeating the argument for the I\ integral yields the existence of some constant 
C a > such that 



lim sup- — — 

^oo log(/?) 



< C a 2 



2S 2 



1 - 5 2 ' 



(236) 
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Since this is true for all 5 > 0, we have 

limsup 

0_>oo \og(/3) 

finishing Step 4. 



limsup ^ = 0, (237) 



Therefore for \a\ < 2b, by Step 4 and Eq. (161 



r I aoo (A:B) 2Jn 

hm sup — - — — < hm sup — -— (238) 

^^oo lOg(^) /3^oo 7Hlog(/3) 



2a 

< 



7T 2 log(2) (min {1, a}) 2 { 2 (mm 
completing the proof of the low-temperature limit. □ 



27 4 1 

in{l,a})'jT^| J ' 



